Novel Theory for Topological Structure of Vortices in a Bose-Einstein condensate 
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By making use of the ^-mapping topological current the- 
ory, a novel expression of V x V in BEC is obtained, which 
reveals the inner topological structure of vortex lines charac- 
terized by Hopf indices and Brouwer degrees. This expression 
is just that formula Landau and Feynman expected to find 
out long time ago. In the case of superconductivity, the de- 
composition theory of (7(1) gauge potential in terms of the 
condensate wave function gives a rigorous proof of London 
assumption, and shows that each vortex line should carry a 
quantized flux. The (/!>-mapping topological current theory of 
V xV also provides a reasonable way to study the bifurcation 
theory of vortex lines in BEC. 
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It is well known that, as semi-phenomenological sce- 
narios of low dimensional BEC continuum, the Gross- 
Pitaevskii (GP) equation and the Ginzburg-Landau (GL) 
equations are of great importance. For neutral super- 
fluid, GP equation is given by M 



IV'IV', (1) 



and the velocity field coming from the current J* = pV^ 
{p — IV'I^) is defined as 
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where i = 1,2,3 denotes the 3-dimensional space coordi- 
nates. For superconductor, the GL equations are known 
as §, 

-^H;ia, --^,)V + aV' + MV^N = o, (3) 
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where the current J' is covariant under U{\) gauge trans- 
formation 



r = epV 



(5) 
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with Ai denoting the external magnetic vector potential. 
In GL theory the velocity take the same form as Eq.(^; 
m and e should be regarded as the effective mass and the 
effective electric charge (especially for Cooper pair, m 
and e should be replaced by 2m and 2e respectively) . In 
all these formulas ip denotes the order parameter, i.e., the 
condensate wave-function, which is a section of complex 
line bundle. 

In theoretical and experimental studies, the curl of V 
is paid much attention to. For a long time, up to now, 
the wave-function is usually expressed in the form 
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then V becomes the gradient of a velocity potential Q{x) 
{V = ^V8), which directly leads to a trivial curl- free 
result: 



V X y = 0. 



(7) 



But nearly half a century ago Onsager and Feynman 
found that this statement must be modified, and Lan- 
dau predicted (5-functions in it, namely, V x V can be 
non-zero at a singular line, the core of a quantized vor- 
tex line 1^. Therefore, it is indispensable to study: what 
is the exact expression for V x V in topology theory? 

In this paper, based on our 0- mapping topological cur- 
rent theory [|j a novel and precise expression for V x is 
obtained, which is just the topological current with inside 
the (5-function of the order parameter. Thus isolated vor- 
tices in BEC, i.e., the topological excitation, can be nat- 
urally created from the zero points of condensate wave- 
function, and be characterized by the quantum numbers: 
Hopf indices and Brouwer degrees of 0-mapping. Using 
the C/(l) gauge potential decomposition theory, the com- 
posed intrinsic electromagnetic gauge potential in terms 
of the wave-function is studied; so a rigorous proof of 
London assumption (Vi = i^Ai) is given, and the essence 
of this relation is revealed. A step further, the 0-mapping 
topological current theory also provides a reasonable way 
to study the spatial bifurcation of the vortex lines, includ- 
ing intersection, splitting and mergence. Being different 
from the others, the bifurcation theory of this paper does 
not need to deal with the concrete form of the wave func- 
tion. At last, it should be pointed out that all the con- 
clusions of this paper do not matter with the concrete 
form of the nonlinear terms in Eqs.(|l|) and (|3|); the non- 
linearity may even be generalized to a form /(I^Al^) [D, 
and all the conclusions are the same. 

As a matter of fact, by means of the (/)-mapping theory 
great progress has been made in studying the topological 
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invariants and the topological structures in many math- 
ematical and physical topics besides here [H-ESl . 



I. (/^-MAPPING TOPOLOGICAL CURRENT 
THEORY FOR V xV 

The basic field of condensate wave- function ip{x) is 
a section of complex line bundle, i.e., a section of 2- 
dimensional real vector bundle on R^: 



(8) 



Following (/>-mapping theory a 2-dimensional unit vector 
is defined as = -[^ (a = 1, 2), where || cj) \\^= (j)°'(j)<' = 
Substituting these formulas into Eq.(|2[), it is easy 
to find out = ^^abn'^din'' , and the curl of V can be 
expressed in terms of n": 



(9) 



Using diu"' = -|- and the Green function 

relation in (f>-space, ^^g^ln||0|| = 27r(5^(0), one can 
directly prove a novel expression for \7xV 
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where = eabdj(j)°'dk(j)^ is the Jacobian vector, 

and 



1 
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is just a simple 2-dimensional case of the A^-dimensional 
(/)-mapping topological current p3{ . This formula includ- 
ing 5^(0) to describe the singularities of (j) is just the pre- 
cise topological expression for V x y that Landau and 
Feynman expected to find out long time ago. Therefore 
an important conclusion is reached: V x = 0, if f 4> 7^ 
0;Vx y ^0, iff$=0. 

The implicit function theory shows that , under the 

regular condition D ^ 0, the general solutions of 

^\x,y,z) = 0, 4,^ix,y,z)^0 (12) 
can be expressed as 

X ^ Xj{s), y ^ yj{s), z ^ Zj{s), {j ^ 1,2, ■ ■ ■, N) (13) 

which represent N isolated singular strings Lj with pa- 
rameter s. These strings are just known as the vortex 
lines. 

In 5-function theory |p5|, one can prove 
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where D(|)s, = (^e^'cmn^lt^), and is the jth 
planer element transversal to Lj with local coordinates 
[u^ , V?). The positive integer (3j is the Hopf index of 0- 
mapping. Meanwhile it can be proved that the direction 
vector of Lj is 



yds I ^, X u ' 



(15) 



Then from Eqs.(p^) and ( [T5| ) we find the important inner 
topological structure of V x F: 
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Vx v>=-y/3. 



where the positive integer (3j is the Hopf index of 0- 
mapping, and r]j is the Brouwer degree, rjj — ±1. And the 

winding number of 4> around Lj is Wj — (3jrjj. Therefore 
the vorticity of vortex line Lj is Fj = Jj, V x • ds = 

■^Wj, where Ej is the jth planer element transversal to 
Lj; and the total vorticity on a surface E should be 



V X y • ds^ — VVF,-. 
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(17) 



We stress that there are no hypothesis in the deduction 
above. Eqs.(0) and (0) are called the differential forms 
of the quantization condition, which cannot be derived 
from the single-valued principle of wave function, and 
are more essential than the integral form (Eq.(|ll" 
For the GL theory, Eqs.(^) and (||) lead to 



A + A^V X B 



-v. 



(18) 



where A is the penetration depth, A^ = In Lon- 

don approximation, p and therefore A, are treated as con- 
stants; hence when noticing i? = VxA, V-i? = and 
Eq.(|lO|), we find a topological equation for B: 

B-\^V^B^—V xV^—6^i^)D{^). (19) 

e ex 

This formula directly leads to 



B-X'S/^B = ^n"^W, j ^S^{x-Xj{s))ds, (20) 



where $o=-^ is the unit flux quantum. We see that in 
simple case Wj = 1, the above equation is just the so- 
called modified London equation ||l6|,|l0| . This expression 
says that, when the condensate wave function ip has no 
zero values, $ ^ 0, i.e., 6'^{$) = 0, and B - A^V^i? = 0, 
which just corresponds to the Meissner state; while in 
the case of mixed state, (/> possesses N isolated zeros, 
(5^(0) 0, thus a type-II superconductor is penetrated by 
an array of N vortices, with each one carrying a quantum 
flux proportional to the winding number Wj . 
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II. DECOMPOSITION OF U{1) GAUGE 
POTENTIAL AND THE VORTEX WITH 
QUANTIZED FLUX 

The decomposition theory of gauge potential in 
SO{N) and SU{N) gauge theories is now playing a more 
and more important role in theoretical studies, because 
it virtually in put s topological ||6|,|l2[| and other important 
informations jl^l to the gauge potential. In the theory of 
superconductivity, ip is a condensate wave function de- 
scribing the charged continuum, so the covariant deriva- 
tive in U{1) gauge theory is introduced to describe the 
interaction between ip and the electromagnetic field: 



Using Eq.(0) we have Si(V') = <i>QS'^{$)D*{!l), which 
gives the topological structure of intrinsic magnetic field 
Bi{ip); Bi{ip) does not matter with the external magnetic 
field. As before, the zero points of i.e., the singular 
vortex lines in superconductivity contribute to intrinsic 
magnetic field as 



N 



B,{^j) ^ <i>oY,wJ —S^{x-x,{s))ds. (26) 



This leads to an important phenomenon that, the mag- 
netic flux coming from the stationary motion of itself 
is quantized 



D,^ = d,^P~i—A,i;, (i = 1,2,3) 
nc 



(21) 
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where Ai is the magnetic gauge potential vector. The 

complex conjugate of AV' is -D*'0* = diip* + i-^A^tp*. and each singular vortex hue Lj carries a magnetic flux 
And the magnetic field tensor is given by 

f,,=d,Aj-djA. (22) 



Multiplying Diip with tp* and D*tjj* with tp respectively, 
we can deduce the decomposition formula for U{1) gauge 
potential: 
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The above expression Ai = Ai{ip) means that the mag- 
netic gauge potential possesses an inner structure in 
terms of charged condensate wave function ip and ip*. 
The inner structure of Ai (ip) with Eq. ( |2^ ) gives a theory 
that in superconductivity how the stationary motion of 
condensate wave function creates an intrinsic magnetic 
field. This is the important physical meaning of decom- 
position of U{1) gauge potential in quantum mechanics. 

Furthermore it has been proved that the covariant 
derivative part [—{tp*Diip — D*%p*-ip)] corresponds to the 
gradient of a phase factor: {diX) Thus this covari- 

ant derivative part contributes nothing to the field tensor 
fij, so it can be ignored, and 



The above decomposition theory of C/(l) gauge poten- 
tial naturally arrives at the conclusion that, in supercon- 
ductivity continuum, the N isolated singular vortices are 
just N isolated topological elementary excitations carry- 
ing with magnetic fluxoid, while their quantum numbers 
are characterized by topological numbers Wj = pijrij. We 
see that the (/)-mapping topological theory in this paper 
is independent of concrete physical models, that gives a 
profound understanding to the nature of the creation of 
the vortex lines and the flux quantization in BEC. 

Comparing Eq.(p^ with Eq.(H) it directly follows a 
simple relation between Vi and Ai 



V, = —A, 
mc 



(28) 
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which is just the London's assumption IQ. We stress 
that, the essence and the significance of this relation are 
not truly realized until now the inner structure of gauge 
potential is revealed and therefore the stationary motion 
of condensate wave function is naturally related to the 
intrinsic magnetic field. 



III. SPATIAL BIFURCATION OF VORTEX LINES 

The Solution ([l3| ) of Eq. ( p^ is based on the condition 
D {(p/x) ^ 0. When it fails, i.e.. 



It should be emphasized that the above U{\) gauge po- 
tential decomposition theory together with (/)-mapping 
theory has been successfully used to study many other 
topological problems in physics |^ . 

In (/i-mapping theory ^^(-0) can be rewritten in terms 
of as Ai{'ip) = ^eatn'^din'' , and fij becomes fij = 
2^eabdin'^djn^ . Therefore the intrinsic magnetic field 
vector from Ai{ip) is expressed as 

Bi{ip) = 7;£i]kfjk = ^Q-^£ijkeabdjn''dkn''. (25) 

Z ZTT 
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at some points (marked as r*) along Lj, the functional 
relationship between coordinate x and z, or y and z is not 
unique in the neighborhood of r*, because the direction 
of the zero line expressed by 

^ = D^ (^]/D^ (^] \r^,^^D' (^] /D^ 
dz \x J \x J ^ dz \x J \x 

(30) 
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is indefinite at r*. Hence this very point r* is called 
a bifurcation point of the two-component vector in 3- 
dimensional space. 

According to the (/)-niapping theory, the Taylor expan- 
sion of the solution of Eq. ([l^) in the neighborhood of f* 
can be generally expressed as p,SO[ : A{x — x*)'^ + 2B{x — 



x*){z - z*) + Ciz - z*Y + ■ 
are constants. This leads to 



0, where A, B and C 



dz 



2B^ + C = 0or C{^f 
dz dx 



2B^+A = Q. 
dx 

(31) 

The solutions of Eq.(|l]) give different branches of the 
zero lines, i.e., the vortex lines at bifurcation points. In 
following four main cases in the branch process are sim- 
ply discussed (the detailed deduction and figures may be 
found in Ref (lO|j|): 

Case I [A ^ Q): For A = 4(5^ - AC) > 0, from 
Eq.(^Tl) we get two different spatial directions at the bi- 
furcation point 



dx 
dz 



11,2 = 



-B±y'B^ - AC 
A ■ 



(32) 



This is the intersection of two vortex lines of different 
directions. 

Case 2{A^0): For A = 4(5^ - AC) = 0, we get only 
one direction at the point 



dx 
dz 



1,2— 7 



(33) 



which includes three sub-cases: (a) Two vortex lines tan- 
gentially intersect; (b) Two vortex lines merge into one 
line; (c) One vortex line splits into two lines. 

Case 3 M = 0, C ^ 0): For A = 4(5^ - AC) > 0, 
from Eq.(pl^) we have 



dz 
dx 



1,2- 



-B ± y/B^ -AC_ 2B 

c " 



(34) 



There are two sub-cases: (a) Three vortex lines merge 
into one line; (b) One vortex line splits into three lines. 
Case A {A = C = 0): Eq.(^ gives respectively 



dx 
dz 



0, 



dz 
dx 



0. 



(35) 



This case shows that two curves normally intersect at the 
bifurcation point, which is similar to case 3. 

It should be noted that, noticing the continuity of topo- 
logical current j from Eq.(|ll|) {dij'^ = 0), at the bifurca- 
tion point the sum of the topological charge of final vortex 
line(s) is required to be equal to that of the initial line(s) 
for a fixed index j: E/ Pj/Vjf = ■ 
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